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Abstract 

The symmetric top is a special case of the general top, and canonical Poisson 
structure on T*SE{3) is the common method of its description. This struc¬ 
ture is invariant under the right action of SO{3), but the Hamiltonian of the 
symmetric top is invariant only under the right action of subgroup that 
corresponds to the rotation around the symmetry axis of the symmetric top. 
So, its Poisson structure was obtained as the reduction T*SE{3)/S^. Next 
we propose the Hamiltonian that describes the wide class of the interaction 
models of symmetric top and axially-symmetric external held. The stability 
of the levitating Orbitron in relative equilibrium was proved. 

Index terms — Poisson reduction, symplectic leaves, 2-form of Kirillov- 
Kostant-Souriau, symmetric top, relative equilibrium, levitating Orbitron, 
Energy-momentum method. 
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Chapter 1 

Reduction of T*SE{3) to the 
Poisson structure for a 
symmetric top 

1.1 Introduction 


The model of the Lagrange top has a long history and it is still the subject of 
the investigations [1]. The widespread use of the group-theoretic methods of 
Hamiltonian mechanics gave the new impulse for the study of this classical 
model [1, 2]. The effectiveness of this model emerged in the study of the 
stability of the magnetic dynamical systems [3, 4], As it will be seen below the 
Hamiltonian reduction from the general asymmetric body to the symmetric 
top leads to the Lie-Poisson structure, embedded in se(3)*. The symplectic 
leaves of this structure are the orbits of the coadjoint representation of SE{3) 
group. 


1.2 Hamiltonian formalism on T*SO{3) 

1.2.1 Representation of the right trivialization for T*SO{3) 

Let’s look at some useful relations for group SO{3) and its cotangent bundle, 
many of which can be found in [5, 2], 
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Group SO{3) is formed by orthogonal, unimodular matrices R i.e. = 
R^^, det(R) = 1. Correspondingly the Lie algebra so(3) is formed by 3 x 3- 
antisymmetric matrices with Lie bracket in the form of matrix commutator. 
Let’s consider a vector space isomorphism": —)■ so{3) is such that [2, p. 
285]: $,ki = -Bikini, = -Isiki^kh where Sm — Levi-Civita symbol. Then 


It/ = ^ X ry; 

^ = I X 77; 

{tv) = 

B^B"^ = ^ 


The scalar product introduced above allows us to state the equivalence of 
the Lie algebra and its dual space so{3)* ~ so{3). Symbol dehnes 
diffeomorphism. As a rule, the diffeomorphisms used below have a simple 
group-theoretical or differential-geometric sense that is explained in the cited 
literature. 

In the representation of the right trivialization which is shown in [1] and 
[5, p. 314] corresponds to the inertial system, we have: 

TSOiS) ~ ^0(3) X so(3), T*S0{3) ~ S0{3) x so(3)L 

Then the right and the left actions of group S'0(3) on T*S0{3) (Cotan¬ 
gent Lift [2, p. 166]) have the form 

f Rg : (R, e T*S0{3) ^ (RB, ^), Be ^0(3); 

[Lg : (R,^) e T*S0{3) (BR,B^B^i) = (BR,Adg-i^) 


1.2.2 Svmplectic and Poisson structures on T*S0(3) 
[5, 1, 6, 7] 

Liouville form on T*S0{3) ~ S0{3) x so(3)* looks like 

( 1 - 2 - 2 - 1 ) 

where dR = ((iR)R“^ is the right-invariant of Maurer-Cartan 1-form. 

Then by using the Maurer-Cartan equation [2, p. 276] for the canonical 
symplectic 2-form we have 

_ _^qT*so{3) ^ ^ ^ (1.2.2.2) 
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^ d/Tli -\- ^TTiSijk^Rjs R ^Rsk 

Any given symplectic structure defines a Poisson structure on the same 
manifold as follows 

{F,G}{z) = n{ip{z),iG{z)) = = -%G, (1.2.2.3) 

where equation = dG is performed for vector held ^g- 

By considering the elements of the matrix R and the components of 
momentum tt as the dynamic variables on T*SO{3), we can obtain the fol¬ 
lowing set of Poisson brackets that completely dehnes a Poisson structure on 
T*SO{3): 

{RijjRkl} 0, ^Tli, Rjk^ EijiRik, -^TTiyTTjy EijiTTi- (1.2.2.4) 

Note that in the inertial system Poisson brackets for the matrix elements 
R are grouped in columns. For example, Poisson brackets for the elements 
of the 3d-column is expressed through the elements of the 3d-column only. 

1.3 Reduction of the Poisson structure for a 
symmetric top 

Poisson structure (1.2.2.4) is invariant under right translations of constant 
matrix B G S'0(3): 

■((RR)jj, (RR)fc;} 0, i^RR'jjn} ^ijl(^RB')lny 

(1,3,1) 

Subgroup G SO(3) is of interest for the symmetric top where = 
{Z G ^0(3) : Zj 3 = 5 * 3 } (it is assumed that the axis of the body symmetry is 
directed along the vector E 3 ) in the body frame). Then (RZ)j 3 = RijZj^ = 
Ris, i.e. 3d-row of the matrix R remains invariant under right translations 
corresponding to subgroup S^. 

Let’s consider projection SO{3) on the sphere (as a set of unit vectors 
ly with = 1 ) 

r ; R H- = RisBi, Ui = Ray (1.3.2) 

where are basis vectors for inertial system. 

This projection generates a map 

r : T*SO{3) 3 (R, tt) ^ (r(R), tt) = (i/, tt) G iPi, (1.3.3) 
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where Ki ~ x so( 3 )* C x so( 3 )* ~ se( 3 )* = (M^(s)so( 3 )*). 

On se( 3 )*, as on any space that is dual to the Lie algebra, there exists a 
canonical Lie-Poisson structure [2, p. 425]. In this case [2, pp. 491,367] this 
structure is determined by the following Poisson brackets 

{o, Vk] = 0, {vTi, Uj} = EijiVu {tT*, TTj} = EijiTTi. (1.3.4) 

Comparing (1.3.4) with (1.2.2.4) shows that a surjective mapping r is 
poissonian. Thus the conditions of Theorem 10.5.1 [2, p. 355] are satisfied, 
and hence. 

Ad ~ T*SO{‘i)/S^ (1.3.5) 

If Hamiltonian H on T*SO{3) is is invariant, then there is a Hamilto¬ 
nian h on Ki that exists in such a way that H = ho f and with trajectories of 
the dynamic system with Hamiltonian H are r-associated with trajectories 
for Hamiltonian h [2, p. 355]. 


1.4 The structure of symplectic leaves for 
the dynamics of a symmetric top 

Proposition 1. Let’s consider the functions on se(3)* 

Ci{u,Tz) = C2 {v>,tz) = u ■ Tz. (1-4.1) 


1. Functions Ci, C 2 are independent in all points se(3)* for each 7 ^ 0. 

2. This means that a common level for functions (l.f.l) when 

Vq ^ Q is a submanifold in se(3)*. 

3. There is such coordinate system (c^, c^, ..., in the neighborhood 

of an any point {i/q, ttq) E se(3)*(iyQ > 0) that in this neighborhood L(iy, 7 r) is 
determined by eguations 


cl = Ci(i/,7r); 
c 2 = C2{u,tz)] 


(1.4.2) 


Remark 1. Coordinates (m1,...,m^) are internal on the submanifold 
determined by the fixed values {C,C). 

□ 

1. Functions are independent in the neighborhood of point {ly, tt) if dif¬ 
ferentials of these functions are independent at this point. 
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We have 


r dCi = Ip-dz/fc + l^dvTfc = 2i^fcdz/fc; 

I dC2 = |r^dz/fc + l^dvTfc = nkduk + z/fedvr^; 

Let’s assume the contrary, there are such a, (3 that adCi + (3(102 = 0, 
then when equating to zero the coefhcients of the differentials dvr^, we get 

> 0)&(/9z/fcd7rfc = 0) —)> /3z/fc = 0,V/c —>■ /3 = 0 

Correspondingly, equating coefficients to zero when the differentials are 
dz/^, we get 

> 0)&(adCi) = 0 —;■ az/fcdz/fc = 0 —)■ = 0,VA: —)■ a = 0 

2 . As it was proved in item 1 the map C* = Ci x 6*2 is regular in points 
(zz, tt) G se{3)* , where zz ^ 0, that is why item 2 is a direct corollary fact of 
Submersion Theorem [11, p. 175]. 

3. This statement is the simple consequence of the Local Onto Theo¬ 
rem [11, p. 175] for mapping C — Co- 


Coadjoint action of group SE{3) on se(3)* looks like ((14.7.10), [2, p. 
492]) 

= (A[zz],a X A[zz] + A[7r]). (1.4.3) 

Proposition 2. Let K = se(3)*\ {(zz, tt) : zz = 0} that is K is an open 
submanifold se(3)* with Lie-Poisson structure (1-3-4) o.nd coadjoint action 
(1-4-3) induced from se{3)*, then 

1. Functions Ci,C 2 are the Casimir functions on K. 

2. is orbit 0(i^Q,,ro) of the coadjoint group presentation SE{3)- 

3. ,ro) the symplectic leaf K, and every symplectic leaf K can be 
represented as L(i^Q ,ro)- 

4- is a Poisson submanifold in K 

□ 

1. Functions (1.4.1) are invariant relative to the coadjoint action (1.4.3), 
so on the basis of proposition 12.6.1 (see [2, p. 421]) we can assume that 
Cl, (72 are the Casimir functions. 

2. For any point (zz, tt) G ^ro) there are such parameters (a. A) that 
(zz, tt) = Ad*^ ^^_i(zzo, ttq). Indeed, zz g and there is a rotation A such 
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that u = A[i/o] and a suitable vector a can be found with the formula 
a = X (tt - A[7ro]). 

3. As it follows from Corollary 14.4.3 and Definition 10.4.3 [2, p. 

477, (i)] the connected component of the orbit of the coadjoint presentation 
is a symplectic leaf. The orbits of group SE{3) are connected. 

4. According to Proposition 1.2 ,ro) is not only immersed subman¬ 
ifold, but also Poisson submanifold ([2, p. 347]) K (and se(3)*). 

■ 

In sections 14.7 [2] and in section 4.4 [ 8 ] the orbits of the coadjoint pre¬ 
sentation of group SE{3) were investigated. Their structure as the sym¬ 
plectic manifolds with 2-form of Kirillov-Kostant-Souriau (KKS) i® 

described in Theorem 4.4.1 (see [ 8 , p. 142]). 

We are interested in orbits with uq ^ 0. Asa consequence of this 

theorem, we obtain for these orbits: 

Proposition 3. Orbit with uq ^ 0 is diffeomorphic to the cotan¬ 
gent bundle of the sphere and the symplectic form on the orbit differs 

from the canonical symplectic form on the cotangent bundle of the sphere in 
so-called magnetic term. 



where p : —)■ is the projection onto cotangent bundle, 

B — 2-form on sphere, v G ’I ^ 


1.5 Hamiltonian formalism on T*SE{3) 


1.5.1 Bundle of orthonormal oriented triads as 


the configuration space for rigid body dynamics 

Let’s consider a frame of reference, associated with the body that is equivalent 
to an orthonormal frame (triad) with the origin x E E^ in the center of mass 
of a rigid body and unit vectors Ei, directed along the principal axes of inertia 
tensor with the same orientation as the inertial system {e)}, (ei x 62 = 63 ). 



That means that the conhguration space for a rigid body coincides with 
the bundle 0 ^{E^) A E^ of the oriented orthonormal 

frames in the Euclidean space E^. 

The elements of 0^{E'^) are 2 ; = (x, {E*}), {Ei x E 2 = E 3 ), then g{z) = 

X e E^. 

On 0~^{E^), as on the principal bundle with structure group SO{3) a 
canonical right action of group S'0(3) is determined: 

re : ^ = (x, {Ei}) G 0+iE^) ^ (x, {B\Ek}), B G ^0(3) (1.5.1.1) 

In addition, there is left action of SE{3) for flat Euclidean space E^. It 
is not canonical in terms of the general theory of the hber bundles. 

Let’s choose a hxed point O G and a hxed triad, forming a Cartesian 
system of reference Zq = (0 ,{ej}), then each point x E E^ is presented by 
a radius-vector x, and each rotation A G SO{3) is presented by matrix Aki 
such that Aci = AkiCk- 

Then the left action of SE{3) for 2; G 0~^{E^) looks like 

l{a,A)Z = {a + Ax,{AEi}) (1.5.1.2) 

This action SE{3) is simply transitive, i.e. any element z G 0^{E^) can 
be obtained by the left action (1.5.1.2) from zq G 0~^{E^) by the only way 

Z 1{x,'R.)Zq, Rik {Rky ('i) (1.5.1.3) 

Mapping 

T : z = (x, {Ei}) {x, R), Rik = {Ek, efl (1.5.1.4) 

is a global map of the bundle 0^{E^), actions (1.5.1.1) and (1.5.1.2) in this 
map are the right and left translations on SE{3) correspondingly. 

TorBoT"^ = RB, °l{a,A)°'^~^ = L(^a,A) (1.5.1.5) 

Diffeomorphism T allowes to consider group SE{3) as a conhguration 
space for the rigid body dynamics that will be assumed further. 
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1.5.2 Poisson and symplectic structures on T*SE{3) in 
the inertial system 

Take into consideration that group SE{3) as a manifold is a direct product, 
we have 

T*{SE{3)) ~ T*(R^ X ^0(3)) ~ T*(R^) x T*{S0{3)) (1.5.2.1) 

Representations (1.5.2.1) are a direct product [9, p. 81-82], where trans¬ 
lational and rotational degrees of freedom are separated from each other in 
the symplectic and Poisson structures. 

(1.5.2.2) 

= dx'- A dpi - ^EijkSRjk A dni + ^TiiEijkSRjs A dRsk 

Poisson brackets that correspond have the form: 

\j^iiRjk\ ^ijlRlk- (1.5.2.3) 

Poisson brackets (1.5.2.3) explain that the basic dynamic variables refer 
to the inertial system. 

Let — Cotangent Lift of the right translation R^ on group SE{3) 
by elements (0, B) G S'0(3), L^Ib) — Cotangent Lift of the left translation 
respectively. 

In the inertial system we have, expanding (1.2.1.2) 

^(h,B) e SE{3), (R,7r)) G T*SE{3); 

Rb ■ (®,R.) (a3,RB); 

■* -^(a,A) : (a^,R-) —^ (a-I-Aa;, AR); (1.5.2.4) 

Rb ■ {{x,p), (R,7r)) ^ ((a;,p), (RB,7r)); 

Ma,A) ■ ((a +Aa;,Ap),(AR,A7r)) 

In order to make a transition for the reference frame that is connected 
with the body, it is necessary, in accordance with the last line in (1.5.2.4) to 
make the following canonical transformation of the dynamic variables. 

r p = R-V; 

<n = R-^7r; (1.5.2.5) 

[ A = 03, R = R 
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where p is a body momentum in the inertial system; 

TT is a self angular momentum of the body in the inertial system; 

P is a body momentum in the body frame; 
n is a self angular momentum of the body in the body frame. 

If we insert expression (1.5.2.5) into (1.5.2.3), we get the following Poisson 
brackets (only nonzero) in the body frame 


Pj} Pp) ^ijkPki 

-(IljjIIjj- -(Ilfc, Pjj} EkjlRil- 


(1.5.2.6) 


Remark 2. We will write the arithmetic vectors that represent the phys¬ 
ical vector components in the inertial system in small bold letters and the 
components of the same vector in the body frame in capital bold letters. For 
example, 






TT = 

(e2,vf) 

, n = 

(^2,7f) 


.(e3,vf))_ 


.(^3,7?)). 


Remark 3. Expressions (1.5.2.2,1.5.2.3) for the symplectic and Poisson 
structures do not correspond to the right trivialization T*SE{3), described 
in [10], the ng'/it trivialization in these expressions affects only the rotational 
degrees of freedom from T*SO{3), and (1.5.2.5) fully corresponds to the left 
trivialization T*SE{3) [10]. 

This is coordinated with understanding of 0~^{E^) as a conhguration 
space for rigid body dynamics and the roles of the right and left actions 
(1.5.1.1,1.5.1.2). 

Remark 4. The map of bundle 0~^{E^) corresponding to the inertial 
system has the advantage that translational and rotational degrees of freedom 
are divided and it is fully realized in the model of a symmetric top. As 
for general top, the kinetic rotation energy is the left-invariant function on 
T*SE{3) and therefore body frame is traditionally used. For general top 
we can combine the advantages of these two maps, if we use transformation 
(1.5.2.7) instead of (1.5.2.5) 


n = R ^tt; 

P = p, X = X, R = R 


(1.5.2.7) 


Using of this map is described in [7]. 
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1.6 Reduction of T*SE{3) to the Poisson struc¬ 
ture for a symmetric top 

As in case of T*S0{3), the Poisson brackets for T*SE{3) are invariant for 
the right translations 

Rb ; ((*,p),(R,7r)) ^ ((a;,p),(RB,7r)), B G SO{3). (1.6.1) 

As we can seen from (1.6.1), variables x,p for translational degrees of 
freedom are not subjected to transformations, they play a passive role in the 
reduction process and they can be omited in further discussion. 

In general case kinetic energy of the rigid body is the left-invariant, but 
not the right-invariant. However, for the symmetric top, the system is in¬ 
variant concerning the right translation from G SO{3), where group 
of body symmetry (it is assumed that the axis of body symmetry is directed 
along vector E^). 

As in Section ^ the conditions of Theorem 10.5.1 [2, p. 355] are satished, 
and therefore, we get 

Pi = T*SE{3)/S^ ~ X KiCP = T*R^ x K. (1.6.2) 

where Pi is the Poisson manifold with such Poisson brackets 
{Xi,Pj} = 6ij, {Ui, Z/fc} = 0, {vTi, Oj} = EijiUi, {vTi, TTj} = EijiTTi. (1.6.3) 

In order to realize a system reduction completely, it is necessary to trans¬ 
form the standard Hamiltonian for a rigid body into inertial system, namely, 
kinetic energy contribution of proper body rotation. This is not difficult to 
perform for the symmetric top, where two momenta of inertia are equal in 
the body frame. 

Using Ii = I 2 = /_L, after a few transformations we obtain 
Tspin{{{x, p), (R, Tt))) = (jy, 77)2 ^ y g 

That means that the Hamiltonian of the symmetric top takes the form 
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after discarding Casimir function (i/,7r)^. Where M is a body 

mass. 

Hamiltonian h depends only on the dynamic variables on Pi. 

The dynamic system is finally reduced to (Pi, h), as according to 
Theorem 10.5.1 [2, p. 355] the dynamic trajectories of the original system 
are projected onto the dynamic trajectories of the reduced system by Poisson 
mapping 

T*SE{?,) 3 ((*,p),(R,7r)) ^ ((aj,p),(r(R),7r)) G Pi. (1.6.6) 

As for the structure of the symplectic of the symplectic leaves of 

the Poisson manifold Pi, taking the results of Section 4 (see Proposition 
3) into consideration, we have 

X (1.6.7) 

1.7 Momentum map for the action of group 

SE{3) 

1.7.1 Momentum map for the action of group SE{3) 
on itself by left translations 

Let’s consider the relations that are important for understanding a momen¬ 
tum map for group SE{3) and for understanding representation of the inertial 
system. 

Let’s give the invariant form for Poisson brackets in the inertial system. 
Let’s assume that v and uj are the constant arithmetic vectors, i.e. their 
components are constant in the inertial system. 

We have 


6ij, 

{ 77 ^, 77 ^} Rjk} ^ijlRlk 

( 1 . 7 . 1 . 1 ) 

Hence 


{x, (p, v)} = v 

(1.7.1.2) 

We also have 


{R, (tt, o;)} = a)R 

(1.7.1.3) 
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Similarly, we have 


{tt, (tt, lj)} = X tt = Qtt 


(1.7.1.4) 


The relations (1.7.1.3) and (1.7.1.4) can also be written in the form using 
physical vectors 

{Ek,{7r,uj)} = u X Ek, u = Uiei (1.7.1.5) 

{tt, (tt, Lj)} = u; X TT (1.7.1. 6 ) 


Expressions (1.7.1.5),(1.7.1.6) show that tt is the angular momentum of 
self-rotation of the body as it generates only rotation of the vectors describing 
self-rotation of the rigid body. 

To generate the rotation of the all physical vectors of the system we should 
take the total angular momentum of the system (i.e. including the orbital 
angular momentum). 

j = 7T + X X p (1.7.1.7) 


Then we get 

'{x, = w X f; 

^ {P, =i^y<F, 

{Ek, = l 3 X Ek, 

.{t 0',^^)} = w X 7f; 


(1.7.1.8) 


From equetions (1.7.1.8) and (1.7.1.2) we find the momentum map that 
corresponds to the action of group SE{3) on itself by left translations. 


Jl : {{x,p),(R,tt)) ^ {p,7r + X X p) (1.7.1.9) 


1.7.2 Momentum map for the action of group SO(3) 
on SE{3) by right translations 

Let’s hnd the momentum map corresponding to the right action of group 
SO{3) on ^E(3). 

Let n be the constant arithmetic vector, that is Q has the constant com¬ 
ponents in the body frame. 

Then 

{R, (n,f2)} = Ril (1.7.2.1) 
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{-Efc, (n, 17)} — w X uj — fliEi (1.7.2.2) 

Despite the similarity with formula (1.7.1.4) the Hamiltonian (11,17) gen¬ 
erates such a rotation that the angular velocity in the rigid body, not in the 
inertial system, is constant. In this case the angular velocity vector is under¬ 
goes rotation: uj = Rikfl^Ci. 

Also we have 

{n, (n,i7)} = -i7xn = n^i2 (i.7.2.3) 

It is in this way that the arithmetic vector 11 should be changed in oder 
for a physical vector vf = vr^efc = HkEk to get 

{7f,(n,17)} = 0 (1.7.2.4) 

Equations (1.7.2.1),(1.7.2.4) show that momentum map corresponding to 
the right action of group SO{3) on SE{3) is 

jR = n (1.7.2.5) 

For a general top the kinetic energy is the left-, but not the right-invariant 
function, so the components of this moment will not conserved. In the case 
of the symmetric top of component Ha is conserved, and in case of reduc¬ 
tion (1.6.2) to the symmetric top this component transfers into the Casimir 
function C 2 in accordance with the general rule (see 12.6.1 [2, pp. 421-422]). 


1.7.3 Momentum map for the action of group SE{3) 
on a Poisson manifold P 

For left translation on T*SE{3) in the inertial system we have 


Lfa,A) ■ ((a + Aa;,Ap),(AR,A7r)) 

(1.7.3.1) 

As 


r : R H- 1 / = RisBi, Ui = Ris 

(1.7.3.2) 

then 


t(BR) 3 (BR)j3 Bj^kRks Bikk'k 

(1.7.3.3) 

that means 


r(BR) = B[i/] 

(1.7.3.4) 
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(1.7.3.5) 


Consequently, we have for the left action SO{3) on P 

/(a,A)((®>P)> (^>^)) = ((a +Aa;,Ap),(Ai/,A7r)) 

Momentum map for action (1.7.3.5) looks like (1.7.1.9) 

Jl : {ix,p), (i^,7r)) 1 -^ (p, TT + £C X p) (1.7.3.6) 


1.7.4 Momentum map on a symplectic leaf ^(j^o,7ro) 

First of all, the Casimir functions C'i,C 2 are invariant with respect to the 
action (1.7.3.6) of the previous section, and, consequently, they conserved 
symplectic leaves, and the generators of this action will be tangential to 
these leaves. 

Therefore we should compare the momentum map on the manifold P with 
its Poisson structure with the momentum map on the symplectic leaf with the 
inner Poisson structure of the leaf that corresponds to the KKS-symplectic 
structure. 

Let ^ = {v,Lo) G se(3). The vector held has the form 

^k{{x,p), (i^,7r)) = ((i; To; X x,u; x p), (cu x u,u} x tt)) (1.7.4.1) 

Momentum map is determined from the relations (see (11.2.3), [2, p. 
374]) 

{F, JK)} = (1.7.4.2) 

Let’s consider the inner Poisson structure on a symplectic leaf (1.6.7). If / 
function on the manifold P, then we can denote / contraction / on subman¬ 
ifold (1.3.2), that is / = Then according to Proposition 10.4.2 

[2, p. 346] {f,g} = {f,g} for example 

(1.7.4.3) 

Using basic dynamic variables x^,pi,P,7ri and Poisson brackets (1.7.4.3) 
as F in the expression (1.7.4.2), we hnd that dynamic variables J{^) satishes 
(1.7.1.8), that means it is moment for a symplectic structure on A(,^q ,^ 0 ) ^md 
the following statement is true. 
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Proposition 4. Let J : P ^ se{3)* is a momentum map, that cor¬ 
responds to the action of group SE{3) on the Poisson manifold P, and 
J{uo,TTo)--^{i^o,Tro) se{3)* is a momentum map, that corresponds to the ac¬ 
tion of group SE{3) on the symplectic /ea/-n-Q) then 
1. Jl ; ((a3,p), (R,7r)) (p,7r + a; X p). 


1.8 The motion equations of a symmetric top 

in the external field 


We have the Poisson brackets 


{Xi,Pj} = 6ij, {z/j, Ok} = 0, {vTi, Oj} = EijiOi, {TTi, TTj] = EijiTTi (1.8.1) 


and Hamiltonian 

tt))) = + + 

Applying / = {/, h} to basic dynamic variables, we get 

p = —V^V(x, ly); 

i > = j^TT X V ] 

TT = {x, u) X U] 


( 1 . 8 . 2 ) 


(1.8.3) 


1.9 Correlation between energy-momentum 
method and Ratiu-Ortega Theorem (see the¬ 
orem 4.8 [12] in the problem of Orbitron 


1.9.1 Necessary condition 

The Energy-momentum method as well as Ratiu-Ortega Theorem formulate 
the conditions in a fixed (reference) point of the suggested relative equilib¬ 
rium. 
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Let’s add the variables {x, p) to the list of variables Proposition 1, item 
3 changing their names to M 5 , Mg,.... This order and name of variables will be 
more convenient in the following calculations without changing anything in 
the essence. It is important that variables c^,c^ jjnumerate/,^ the symplectic 
layers ,^ 0)5 while variables ui, U 2 , ■ ■ ■ are the internal variables on the leaf. 

So, let Zq € A(i^q TTo) C P, (^'o > 0). Then there are 2 coordinate systems 
in the vicinity of this point: the 1 st one is global P = ((i/, tt), {x,p)) and 
the 2 nd, generally speaking, is a local (c^, c^, u^, u^, u^, ...). 

We assume that the augmented Hamiltonian on the symplectic leaf is 
the restriction of the augmented Hamiltonian on P (this is done for the 
Orbitron and in a more general case, for example, for polynomial according 
to (i/, 7 r)) Hamiltonians). 

So, the 1st condition in the Energy-momentum method is 


= 0 

( 1 .9.1.1) 

For any function / in the vicinity zq we have 


d/ = dc/ duf 

(1.9.1.2) 

Condition (1.9.1.1) can be written down as 


d«/if = 0 

“ |2o 

(1.9.1.3) 

It is also obvious that 

f duCi = 0 ; 

|d„C '2 = 0 ; 

For any constants A^, A^ is performed 

(1.9.1.4) 

d„(h« + AiCi + A2C2)po = 0 

(1.9.1.5) 


In order to fulhll condition (1.9.1.5) not only for partial but also for the 
total differential, A^,A^ must have specihc dehnite values that satishes the 
equations in point zq, 


X^d.iCi + X^dciC2 = 
X^d,2Ci + X^d^2C2 = 
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Due the special choice of the coordinate system (c^, c^, v}, u^, u^, ...) 

we have 


= -9^2 hS; 

(1.9.1.7) 

If A^,A^ satisfy (1.9.1.6),(1.9.1.7) then 


d'li = 0 

(1.9.1.8) 

where 


= h^ + X^Ci + A^Gs 

(1.9.1.9) 


Condition (1.9.1.8) has the form of 1st condition in Ratiu-Ortega Theo¬ 
rem. 


1.9.2 Sufficient conditions 


Let’s consider the 2nd condition of G^-stability for Energy-momentum method 
and Ratiu-Ortega Theorem. 

First of all, let’s consider Hessian transformation in transition from one 
local coordinate system (y^) into another one {z^). 

Q2f ^ Q2J Qyk Qyl QJ Q2yk ^ l O O l ^ 

dz^dz^ dy^dy^ dz^ dz^ dy^ dz^dz^ 


Hence 


d'^f df 

dlfivi, m) = rj2), vi, m e T,,P 

If the coordinate system y^ means (c^, c^, v?, u^, ...) then 




V c^v 

dc^dc^ 


dc^du^ ) 

(9^ f 

j_ i _V 7;"V ij^ 


(1.9.2.2) 


(1.9.2.3) 
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iVuV2) 


dy>^ 

Now let’s assume that each of the vectors r}i,r ]2 G ,^ 0 )) that means 

. d 


7] = rj 


du’^ 


(V^c^ = 0)&(VX = ^“) 


then 




dy^ 

dz/(hi,^ 2 ) = (1^/(771,772) + (V/,d2y(77i, 772 )) 

Let’s apply formula (1.9.2.6) to function at the point zq 

d^h«’^(77i, 772) = d^h«’^(77i, 772) + d^77(77i, 772 )) 


(1.9.2.4) 

(1.9.2.5) 

(1.9.2.6) 

(1.9.2.7) 


Because of the conditions (1.9.1.8) 2nd term in (1.9.2.7) is equal to 0. 
Therefore we have 


d^h^’^(77i, 772) = d^(/7^ + A^Ci + A^C'2)(??i, ^ 2 ) (1.9.2.8) 

since d^Ci = d^C '2 = 0 then 

d^h«’^(77i, 772 ) = dlh^{r]i, V 2 ) = d^/^( 771 ,772). (1.9.2.9) 

Thus, the positive positive definiteness of the quadratic form in Energy- 
momentum method (the right side of (1.9.2.9) is an equivalent to the positive 
dehniteness of the corresponding quadratic form in Ratiu-Ortega Theorem 
(left side of (1.9.2.9)) on the vectors that are tangent to a symplectic leaf. 

We also note that the requirements set by the Energy-momentum method 
on the subspace W in a case studied coincide with the corresponding Ratiu- 
Ortega Theorem. 
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Chapter 2 

Hamiltonian dynamics of a 
symmetric top in external fields 
having axial symmetry. 
Levitating Orbitron 

2.1 Introduction 

We propose a Hamiltonian that describes a wide class of the models of a 
symmetric top that interacts with external field having axial symmetry. The 
stability of relative equilibrium is considered. It investigates on the base on 
Ratiu-Ortega Theorem. 

We found the necessary and sufficient conditions of the dynamic equilib¬ 
rium. For the first time in general form the final analytical expressions for 
sufficient conditions were deduced. 

The general formulas are tested for two important physical models: the 
model of the generalized Orbitron (see [10] - without gravity force) and the 
model of levitating generalized Orbitron - new task that generalize the results 
of the article [13]. 

Thanks to the new approach became possible to formulate the clear fully 
analytical proof that levitating of the Orbitron is possible. 

The whole study is conducted in the inertial reference system, that is 
justified in the Ch.l of this article. Where the equivalence of the algorithms 
of energy-momentum method and Ratiu-Ortega Theorem is also proved. 
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2.2 The motion equations of a symmetric top 
in an external field 

We have the Poisson brackets 

{Xi,Pj} = 6ij, {Ui, z/fc} = 0, {vTi, l/j} = EijlUl, {TTiyTlj} = EijlTTi. (2.2.1) 

h{{{x,u),{p,%))) = + + 

Applying / = {/, h} to the base dynamical variables, we get 

p = —V^V^x, z7); 
z7 = -^TT X z7; 
vf = {x, z7) X P; 


( 2 . 2 . 2 ) 


(2.2.3) 


2.3 Axial symmetry condition 


By virtue of the fact that kinetic energy S'0(3)-symmetric then the condition 
of axial symmetry of the Hamiltonian can be reduced to the condition of axial 
symmetry of the potential energy. 


^dV _ ^dV 
^ dx'^ ^ dx^ 
that is also can be written in the form 



+ z/^ 


dV 

(9z/2 


— u 


dV 


0 


(2.3.1) 


fx X (V'lx)^ + Pi X (W)i = 0. (2.3.I0) 

In what follows we will show that (2.3.1a) is identically fulhlled on the 
orbit of the relative equilibrium, so it bring nothing new in the study of 
stability. 


2.4 The motion integrals and augmented Hamil¬ 
tonian 

The motion integrals and augmented Hamiltonian form remain the same as 
for the Orbitron, that is 

= h — ojJ^ + AiC*i + X 2 C 2 (2.4.1) 
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where 


{Ci{u,Ti) = z72; 

= (2.4.2) 

[ J 3 (((f, z7), (p,^))) = VTg + XiPa -X 2 Vl- 


2.5 


Necessary condition of the relative equi¬ 
librium 


= I -p 


ooe^ X X 


dp + ( 


ues + X2^ ) ■ 


+ {V^V + uje3xp)-dx + {V^V + 2Aiz7 + Asvf) • diJ 
we getting the necessary condition of the relative eqnilibrium 


( 2 . 6 , 1 ) 


'p = Muj{e3 X if); 

V^V + cues X p = 0; 

TT = I±UJe3 - A2/±zA; 

+ 2Aii7+ A2vf = 0; 


(2.5.2) 


( VW = Mu;^x±] 

I A 2 = oji ^3 - t^C 2 ] (2.5.2a) 

[ V'^14 + Az7+ 1 x 00 X 263 = 0; 

where 

A = 2Ai - Xpx (2.5.3) 

It follows 

r V^I/ = Mu^xx] 

I {V'^V)x = -Xux; (2.5.4) 

[ujC 2 = 0 j{u, ff) = + (A + Ix^‘^)l'z 

Here, the symbol ux is a component of the vector u that is orthogonal to 
the symmetry axis z. 

Remark 1. Vector z7 and Lagrange multiplier oo can be found from hrst 
two lines of (2.5.2) or that is the same from the hrst line of (2.5.2a). That is 
distinct from Orbitron, where we postulate that z7 direct in the line of 2 ; axis. 

Remark 2. Formulas (2.5.3) not only determined z7 and Lagrange mul¬ 
tiplier cj, Ai,A 2 , but most likely applies restrictions on the function V in 
addition to it axial symmetry. 
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2.6 Supporting point and allowable varia¬ 
tions 


'xq = roCi; 

^ Po = Mujroe2] 

^0 = i^O, 1^0 = 1 
^TTo = 6*21^0 + 

As it was already mentioned, vector z/q must be found from the necessary 
conditions of stability. 

Allowable variations, vector v = {6x, 6p, ( 5 z/, 671) in the supporting point 
fit the conditions of 


(2,6,1) 

5 


SvC2 = TToidui + uoidTii = 0 ; ( 2 . 6 . 2 ) 

( 5 ^ J3 = 6113 + po6xi + ro6p2 = 0. 

Adding condition of transversality to these conditions that in contrast 
to the Orbitron task we will take in more simpler form that simplihes the 
further computations. 

So, the complete set of the conditions can be written in the following form 

6773 = —— (i^o±, (Stt — 

^ 0 . \ V ‘'Oz (2.6.3) 

Sp2 = ——Sxi — 

Jx2 = 0; 


Remark 3 . As it follows from second line of ( 2 . 6 . 1 ) that for the vector 
tangent to the orbit the equation 0:02 = 0 could not be fulfilled in case of 
conditions of ro 7^ 0 and a; 7^ 0 that are always supposed. 


2.7 Coordinate system on the plane that is 
fit to analysis of restraints of the varia¬ 
tions 

It is convenient to provide analysis of the restraints of the variations in the 
coordinate system that is fit to the vector z 7 o±. 
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Remark 4. If this vector is 0 (as in the Orbitron) then the new coordinate 
system is coincides with original Cartesian system on the plane. 

Let’s consider the coordinates system on the plane is dehned by the vector 
of a. 

Let on the plane we dehned the coordinate system {ei, 62 } and vector a. 

Let’s direct the basis vector Ei in the new coordinate system {Ei, i? 2 } in 
the line of vector a. In order that new coordinate system has been properly 
oriented {Ei x E 2 = 63 ), the following equations must be satished 



(2.7.1) 


Allow the matrix 


^i""^ 5 Ei OCki^kj 0.ki L/j) 

_|a| |a| _ 


(2.7.2) 


Next 


O^ki^i (Cfc) ^lEi) (Cfc) x') Xk 


i.e. 




(2.7.2a) 



(2.7.26) 


Correspondingly 


Xk = (a ^)kiXi, X = a ^x 


(2.7.3) 


where 



(2.7.3a) 


. |a| |a|. 

In our case, the vector dehned on the plane is so 


-I- I'm - 1^02 

|ih±| .^02 1^01 


, EA = aBAeB, CUBA = {&b-,Ea) (2.7.4) 



(2.7.5) 
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2.8 Restraints of the variations in the new 
coordinate system 


We shall use the reference coordinate system for the variations 6x, 6p and 
the basis {Ei, E 2 ,e 3 } for the variations 6u,67r, and denote these variations 
in the new coordinate system by the capital letters 6N, 511. 

Then restraints (2.6.3) have the form 


5 1/3 
dTTs 

Sp2 

5X2 


I^Q-hl 

Wo±\ 

l^Oz 


5Ni; 


5ni - 




-Mcudxi + ^(6Ui- ; 

^ roi^Oz \ ^ I'Oz ’ 

0 ; 


Allow the variable 5n( 


I 5n; = 5ni - 

\5ni = 5n; + ^/±a;5Ni; 

Then 

(S 123 = -^5Ni; 

Stts = -^ 5 m; 

6p2 = -Mudxi + ^ 5 n;; 

^5x2 = 0 ; 


( 2 . 8 . 1 ) 


( 2 . 8 . 2 ) 


(2.8.1a) 


2.9 The original and reduced quadratic forms 


The second variation of the augmented Hamiltonian has the form (hereinafter 
indices i, j, k = 1..3, A,B,C= 1..2) 


+ ^57f^ + y-'^TTg 

I± Ix 

+2A2(57f_L, Su±) + 2A257r35z/3 
+2Ai5i7| + 2 Ai5z/3 


(2.9.1) 
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-2u{6xi6p2 — 6x2Spi) + 


d^V 

dx'^dx^ 


6x^6x^ 


n . i. A o 3 

+2 ^ . 6 x + 2 ^ ^ „ 6 x 


dx^du^ 


dx'-di'^ 


dv^dv^ 

By using restraints we get 




du^du^ 


dvl 


l,Spi = MuiHxl - 2A^_Sx,m\ + 


-5tt‘^ = 


1 1 


rovoz 

J_ 

I" 0 , 


^^o’^Oz 


5nf; 


6Uf + 2^a;(5n;(5Ni + f 


(57?, (5i7) = ^(5n;(5Ni + 6 U 26 N 2 + ^J±uj 6 Nl, 


5z72 _ ^^^2 ^ ^^2 

^— 2 uj{ 6 xi 6 p 2 — 6 X 2 SP 1 ) = 2Mu‘^6x\ — 2uj^6xi5Yi'^ 
Hereinafter we use the following notation 


«lMI = dW(e-.B^); 

35g5jj = d'V(B*BB); 

aV ^ Q_ (av\ . 

dNAdu3 ^Ea I 9id3 / ’ 


(2.9.2) 


(2.9.3) 


besides that in the expressions of this type we always considered that the 
vectors are tangent to the variables space of x, and Ea to the variables 
space of z7. 

Substituting (2.9.2) into (2.9.1) we hnd the reduced quadratic form 


+ j^Spl (2.9.4) 

+Muj^6xl - 2j-^5xi6U[ + 

roi^oz Mrlul 

+ 2 —a;(5n;(5Ni + + Xul 

fn ^oz ^oz j_i_ 

+2A2 I ^5nUNi + (5n2(5N2 + —/xa;(5N? ) 

V^Oz ^Oz ) 

+2Ai (Xs^l + d^lj 
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+2Muj^6xl - 2j-^6xi6U[ 
ro^oz 

, 2 n r r r 2 

+ 2 dxidx3 + -^dx^ 

OX^ OXiOX^ ox^ 


+2 


1 

d^V 1 d^V 


Wo±\ d‘^V 


i^oxl d'^V 


f, If, - 2 '"'^' , Sx^SN, 

Uqz OX^OIX^ Vqz OX'^OV^ 

-B»:* •«•«"■*&« 


M ^n; 


|i^±| d'^V 


-^7j^7rT'^Ni<5N2 

vqz aN 2 C/i^-^ 


z72^ av 


AN? 


2.10 The method of successive elimination 
of isolated squares 

To derive the conditions of positive definiteness of the quadratic form 

Q = Q{xi,...,Xn) (2.10.1) 

we use the method of successive elimination of isolated squares. 

Let’s represent Q in the form 

Q{xi, ...,Xn)= Ax\ + 2B{x2, ..., Xn)xi + Q'{x 2 , • • •, x„), (2.10.2) 

where B linear function of its variables, and Q' is quadratic function of the 
variables but without xi. 

For the positive dehniteness of Q is necessary that A > 0 then Q positive 
dehniteness is equivalent to the positive dehniteness of the quadratic form 
Qi, but with smaller number of variables. 

Qi{x2,. . . ,Xn) = -^5^(x2,. . . ,X„) + Q'(x2 , . . . ,X„), (2.10.3) 

Consistently applying this procedure, we can hud all the conditions for pos¬ 
itive dehniteness of the original quadratic form, that is means all A > 0. 
Remark 5. The order of elimination of the variables is arbitrary. 
















2.11 Sufficient conditions of stability 

Applying the method of successive elimination of isolated squares in the 
following order 6 p 3 , 6 pi, 5112,511']^, 5N2, 5Ni, result 

' X + (PV{E2,E2) > 0 ; 

A + d'^V {uq , {j^ozoo + A2)^ + - 

A>0, ^>0, AC-B^>0; 


d^ViE2,i!X f 
A+d 2 V(£ 2 ,E 2 ) 


> 0 ; 


where 


A — hdoj 


I d'^V 

AMrl - ^ d'^V (,9^* 


I±i^± + Mrl dxl X + 


92 V 

aN? 


( 2 . 11 . 1 ) 

(2.11.2a) 


r, l±Wo±\po ( ,, I I N I a2y(E fVT'i 

+ A2j + d lei, z/q j 


aN;^ 


d^V(E2,y^ ? 


X + dW (z7([, zA([) + _j_ y^) + i^j_I±uj‘^ - a^v 

8N: 


B = 


I±i^0±+MrS 

d^V 

dxidx^ 


( d'^V A ( d'^V 
I 9xi9N2 j I 9x39N2 


A + 


92 V 

9N? 


( 2 . 11 . 26 ) 


) I±Wo±\po 

‘I±^o±+Mrl 




^ {eqzU! + A 2 ) + d^lA(ei, Eq) — 


dx^dN2 


A+ 


a 2 y 

aN 2 


A + d^l/ {uj,uj) + i^oz^ + X 2 ) + - 


2 _ d 2 V(£; 2 ,^T )2 


X dV(e3,z7([) 




A + 


92 y 


C = 


( d'^V A^ 
d’^V l9x39N2j 


dx\ 


A + 


92 y 

9Nf 


A+ 


d'^V 

dlA'i 


(2.11.2c) 


d2^(e-'3. 


4^d2y(E2,^7T) 


AH" 


a'-^v 


X + dW {uq , 1 X 0 ) + {j^OzOO + A2)^ + z^_Ld±<^^-° ^ 


I±u^^+Mr^ 


A+ 


d'^V 
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We use the following notation. 
Allow the vector 


Uq = - \uo±\e3, (z7o, ) = 0 


'T\ 


(2.11.3) 


Then 


d^V{E2,iyo ) = uq 

) = J^Oz 

) = t'o. 


where 


d^V{ei, 

d^V{e3,u, 


d^V 
^9NiaN2 

a^v 


a^v 


r? I a^v 
^o-L I aN2au^ ’ 


d^V{ua , y, 


iT 


\ ^01. \ Qx^du'^ ’ 
Irr I a'^V 

0 J — ‘'[)Z dx^ai<ii \^0±\gx3du3^ 
■T\ _ ,,2 a^v ouv I,, a^v 


(2.11.4) 


■2 a^v 


'0 ’ ‘'0 


•\ — ,,2 a^v _ 91,7 I,, a^v I ,72 ^ 

J ^OZ^^ + y0± Q^\ 


a^v 

axiSNi 

av 

aa:iaN 2 

a^v 

aNiai^® 

_^v 

aN^al^ 

aW ^ 

a^v 

aNiaN2 
av _ 
Mi 


aW 


i^Ol 


+ 


aW 


axiai/i |i7o±| 
a^ V UQ2 

dxxaux \i'oi-\ 
a'^V izni _|_ 
au\au3 i^oxi 
a^y UQ2 


^02 


axidu2 l^oxl ’ 

I a^y UQi . 

axidu2 li'oxi ’ 

a^y 1202 . 

di'2du3 \u()^ \ ’ 
j a^y 1^01 


duidu3 |i7oxl av2av3 |i7oxl 


lAxP 


aW ,,2 


du^'-'Ol 


+ 2 -: 


1 


I Ax I 


lAxP 
1 f a^v 

du'i 


a^y u 7/ -U a2y,,2 \ . 
MrM;^owo 2 + -^^02) > 

( d^y,, ,, _L d^v,. , a 2 y / 2 

(^“Mf^oWo2 + ^i"0W02 + a;;^Woi “ 

aUI^^oiyo2 + -^^01J , 


(2.11.5) 


02 . 


^02 2 - 


2.12 The case of a dipole in a magnetic field 

Let’s study the possibility of the levitation of a dipole in an axially symmetric 
magnetic held. 

The potential energy in this case has the form 

V = —/i(y, B) + Mgz = —jjiEBi{x) + Mgz (2.12.1) 

where B is an axially symmetric magnetic held. Next, everywhere we assume 
that p > 0 . 

From (2.12.1) follows 

V^l/ = -f^B ( 2 . 12 . 2 ) 
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Hence 


{ A 2 — OJT^z — 7 ~C' 2 ; 

a;C*2 = + (A + (2.12.3) 

From (3) it follows that Fj_ has the same direction as the B± (with the 
sign). As it will be shown below, this condition is coherent with the equation 

= Mu‘^x± (2.12.4) 

for the potential energy of the dipole (2.12.1). 

Therefore, the supporting point in this case can be selected in form 


'xo = roCi; 

Po = Mujroe2] 

i^Q = l^zes + IXrCi, ^1 + 1x1 = 1 , Vr = 

^vfo = I±uje3 - A2/±Fo; 


(2.12.5) 


i.e. vector Fq lies in the plane span(xo,e 3 ). 

In the dipole case the signihcant simplihcation will be 


dlV = 0 


( 2 . 12 . 6 ) 


2.12.1 The relations for axially symmetric magnetic 
field 

It is assumed that the sources of magnetic held are localized, and motion of 
the magnetized body occurs in an area free of sources. 

Then in the movement area the equations of the magnetostatics without 
sources are fulhlled. For the axially symmetric magnetic held in particular 
that means in the cylindrical coordinate system there are only two compo¬ 
nents of the held Br, B^ that depend on 2 variables — r and 2 . 

In given case we have the following equations of magnetostatics 


Br^z Bz^r 0 , 

Bz,z + -Br,r + “-Br = 0; 
Bz,zz T Bzrr T zBz.r 0 


( 2 . 12 . 1 . 1 ) 
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For the Jacobian we have 

' B A, c = Be, A = ^5 AC + {Br,r - f) 

B3,C = 

-83,3 = Bz,z 

For the Hessian of the magnetic field we have expressions 


( 2 . 12 . 1 . 2 ) 


_ _ D 

^z,rz ^3 


Ba,cd 

_1 I 2 Br \ / xaScd+^c^ad+o^d^ac _ 4:Xaxcxd \ . 

r ' r ) \ r / ’ 

-S 3 ,CD = \Bz,r5cD + {Bz,rr “ ^-Sz,r) 

B 3 ,C 3 = ^{Bz,zr); 

-S3 33 Bz ZZ1 


(2.12.1.3) 


(2.12.1.2a) 


In the supporting point where x = (ro, 0, 0)^ for the components of the 
Jacobian we have 

-Si,i Bj-ri 

51 .2 = -S 2 ,i = 0 ; 

52.2 ~ \Br] 

Bz,\ Bz,ri 

53.2 = 0 ; 

, S 3_3 = Bz,z 

For the components of the Hessian in the supporting point we have 

-Si,11 = —Bz,zr + \ {Bz,z + f-Sr) = —Bz,zr ~ ^ {Br,r ~ ^-Sr) J 

51 .12 = -Si,21 = 0; 

51.22 = ^ {Bz,z + f-Sr) = {Br,r “ r^r) ] 

-S2,ll = -Si,12 = 0; 

52.12 = -S2,21 = -Si,22 = “ {Bz,z + f-Sr) = (-Sr-,r ~ ! 

52.22 = 0; 

-S3,11 Bz,rri 

53.12 = 0; 

53.13 Bz,rz) 

5 3 .23 = 0; 

, S3,33 = Bz,zz'l 

(2.12.1.3a) 
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2.12.2 The equations of equilibrium for the magnetic 
dipole 

From the relation (2.12.1) we have 


Ks = + Mg = + Mg] 


c “ Bd^c — - 83 , 0 ; 

dV 
dx^ 


( 2 . 12 . 2 . 1 ) 


Taking in to account the third equation of (2.12.3) the properties of the 
axially-symmetric magnetic held we obtain 


Ur = jBr] 


( 2 . 12 . 2 . 2 ) 


and equations of equilibrium take the form 


+ u'^Br^z = fg; 

U^Br,z + l^'’Br,r = -fuj'^T] 

+ Z/2 = 1 


(2.12.2.3) 


2.12.3 The sufficient conditions of stability for the 
magnetic dipole in an axially symmetric mag¬ 
netic field 

In the case of the magnetic dipole in the supporting point we have the fol¬ 
lowing relations 

axi^Ni Qxidui ’ 

9^V - siffnfi/ )M^- 

9^V - simlz/ i ■ 

9^v = i ■ 

dN2dl'^ ^ ^ ’ 

d^v _ d^v. 

< Wf aiif’ (2.12.3.1) 

av ^ d^v . 

9Ni9N2 dvidv2 ’ 

92 V ^ d'^V . 

9N| ’ 

d^y(ei, ) = sign(z/^) 

^d^y ( 63 , ) = sign(z/^) f - Vr 


d^V \ . 
dx^du^ j ’ 

9V ^ . 
dx^du^J ’ 
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dxidui 

d^V 

dx 181/2 

a'^v 

8 x 381/1 

av 

8x38^2 

a^y 

axiavi 
d‘^V d^V 

^ dx^dui ^ ^ dx^du'^ 


= 0 ; 

fj/Bp^^i 

= 0 ; 

Qx^du^ f^i^Bp^pUz Bp^z^p^ ^ 

t^i^Bp^z^z Bz^z^r) •> 


^ + ^Br) ] 

fi{Br,zl^r - Br,r^^) = Mg+ 
fJ/{B^^^Ur - Br,zl^^) = Mui^r - 


(2.12.3.2) 


(2.12.3.3) 


(2.12.3.4) 


Then the stability conditions of the relative equilibrium in 2.11 can be 
reduced to the form 


A>0, ^>0, AC-B^>0; 


(2.12.3.5) 


A = Moj 


^3Mrl-I 

Jxy 2 + Mr 2 


jJiU^Bz^rr + 


/iU^— I Bz^z H —Bj 


'i!ii+Mrl + ^ 2 ) + {Mg + ^Vz)) 


(2.12.3.6a) 


A + 




^ {uzU + As)^ + y2J_La;2 


'/x 


7]M+Mlg 

B = -llP^Bz^rz - 
M+Mrg + ^ 2 ) + {Mg + 


(2.12.3.65) 


A + r {i^zOJ + \2f + y2/_La;2 




Moj^r — 


fiBr 


{Muj^r - i^Ur)' 


^ + iJt+Mr^ + ^ 2 )^ + y2j^a;2 


(2.12.3.6c) 


where values of the variables i/r,i^z,^ can be found from (2.12.2.3) and 
A 2 ,A, 7 fo from (2.12.3). 
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2.13 The Orbitron 


The most simple example is the Orbitron. 

As the field of the Orbitron (we mean generalized Orbitron see in [10]) 
we assume the field axially-symmetric about an axis 2 ; and mirror symmetric 
respect to the plane vertical to axis 2 ;. 

Let’s show that there are stable relative equilibrium that are located in 
the plane of 2 (equatorial plane). 

In this case without gravity force 

V{x,u) = —/ii?i(T)z/* = —jjL{B{x)^V) (2.13.1) 

In view of the mirror symmetry the supporting point is given by 


and, therefore, 


Xq = roCi; 

Po = Mur 062] 

Uq = 0-63, a = ±1; 

TTo = CTC'2e3 = 7r3e3 = TToCs; 


(2.13.2) 


= ^( 


(2.13.3) 

= a^iBz + Ci;(7ro — I±u) 

'Bz,z 

= 0; 


^ z^r 


(2.13.4) 

B = 




Now we consider the sufficient conditions. 

By using all previous simplifications of (2.12.3.5, 2.12.3.6a-6c) we obtain 


A = + cuTTo — I±_u'^ > 0; 

5 = 0; 


(2.13.5) 


A = 3Mu^ — /laBz^rr > 0 


C = —a^Bz 


[Mu'^r 


2^^2 


> 0 


(2.13.6a) 

(2.13.6c) 
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From (2.12.3.6c) follows —cBz^zz > 0 then the conditions (2.12.3.6a,2.12.3.6c) 
can be written in the form 


-o-Bz,zz > 0 ; 

■O' (S^Bz^r + Bz^rr) = “O' [—Bz^zz + 2^i?2,r) > 0; 


(2.13.7) 


uttq > -ajiBz + ; 

V — (J^Z,ZZ 

These conditions are equivalent to the conclusions of the article [10]. 


2.13.1 Dipoletron 

In standard Orbitron [3] as the sources of the magnetic held the magnetic 
poles were used. The results shown above demonstrates that other kinds of 
the held sources are also allow the stable conhnement. 

A special interest are the sources of the held in the form of magnetic 
dipoles, because of recognized model of a held. Furthermore, it is known 
that stable dynamic conhgurations not realizable in the system of magnetic 
dipoles that interact by magnetic forces only the (so-called ’’problem ^ [14]). 
This circumstance adds the interest to such choice of the model held. 

Let’s consider a system with two magnetic held sources that are the mag¬ 
netic dipoles located on the axis z at points z = ±h and with equal magnetic 
moments oriented along the same axis . 2 . Obviously, the held of this system 
is axially symmetric about axis z and mirror-symmetric with respect to the 
plane z = 0. 

Remark 6. In our system, we have outside forces that keep the dipoles- 
sources in a predetermined position. 

The held of the magnetic dipole is well known 


47r i R^ 



(2.13.1.1) 


where rh vector of the magnetic moment, and R is the radius-vector from 
the dipole to the point of held observation. 

For the dipoles located on the axis 2 ; in the points ±/i in the components 
of Cartesian system we have 


D±_ 2{x3^h)'^-{xl+x‘^) 

=F 


1 , 2 ; 

D± = xl + xl + {X3 =F h)2; 


(2.13.1.2) 
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where q = |^|m| — ’’magnetic charge” equivalent to a pole of a magnetic 
dipole. 

By using (2.13.1.2) and take a derivative of with respect to r at 2 ; = 0 
we can express all quantities of interest. 

B,,r = Br,z = r/3,^ = -6qro{-rl + 

' 5 ,,,, = 6q{34 - 2Arlh^ + 

- + ^ = -6g(r4 - 23Tlh^ + 

where Dq = + h?. 

Then from the first two conditions (2.13.7) we obtain the geometric con¬ 
ditions for the Dipoletron system 


3(a)4_24(m)2 + 8<0; 

(m)4_28(a)2 + i6 > 0; 


or 

2^1^^^^^ < ^ < (2.13.1.5) 

with fulhlling of geometric conditions (2.13.1.5) the third condition (2.13.7) 
can alway be satisfied. 


2.14 Levitation of the Orbitron 

Let’s explore the possibility of the levitation of a dipole in an axially sym¬ 
metric magnetic held. 

As it was shown above, the held of the Orbitron can provide the equi¬ 
libration of the centrifugal force with stability. However, the studies show 
that this held is ill-suited for the equilibration of gravity force. Therefore, 
add the held that linearly depends of coordinates in the system. 

Let 

B{x) = B^{x) + B^{x) (2.14.1) 

where B^ — the magnetic held is linearly dependent on the coordinates, and 
B^ — the held that mirror-symmetric with respect to the plane 2 ; = 0 (of 
course, both helds axially-symmetric). 
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It is expected that the held for the most part is intended for the 
compensation of the gravity force, and held for the compensation of the 
centrifugal force. 

Let’s consider the relative equilibria that are spatially observed in 2 ; = 0. 
The main properties of these helds are 


r bx =^ = 0 ^ (50 = o)&(i?o, = o)i,=o 

\ = 0)U=o; (2.14.2) 

[d'^B^ = 0 ; 


We have (see [15]) 


Then 


B^ — Bq + B'z] 

= -\B'r- 

(2.14.3) 

B^ = B'- 

= B^^, = 0; 

JDL _ 1 TDf. 

(2.14.4) 

[ S... = B'; 

Br,. = B° ; 

Id _ 1 TDf . 

(2.14.5) 


2.14.1 Necessary conditions of the relative equilib¬ 
rium 

Equation of equality of the forces can be written as 


where 


{ Z/^ + /3z/^ = k; 

+ vl = I 



B° 

r,z . 
’ 

Mg . 
uj'^r . 


(2.14.1.1) 


(2.14.1.2) 


38 


From here 


(2.14.1.3) 


From the equation 

Vr = ^Br 
A 

and (2.14.3) we have 


Ur 


Mgr 

2kX 


sign(z/^) 


—sign(K) 


(2.14.1.4) 

(2.14.1.5) 


(2.14.1.6) 


Thus the sign of Ur should be opposite to the sign of k. 

Therefore, in the expression (2.14.1.4) 1-st member and must be negative, 
then the sign of the f3 must be the negative. 


/3 < 0 


(2.14.1.7) 


Furthermore, since vector z7is completely determined of the system (2.14.1.1), 
then the value of A is also dehned by (2.14.1.6). 

Therefore, the equation 

A =/i—+ ^ -/xa;2 > 0; (2.14.1.8) 

should not be regarded as an equation for A, but rather as a limitation for 
Bz and vr. 

From (2.14.1.3) and (2.14.1.4) we have 


/3 1-p/j yi-I-_ ^2 

2(l+/32) I + /32 K 


(2.14.1.4a) 


If |fi;| = 1 then this value with minus sign before the square root becomes 
negative. 

Obviously (because of /5 < 0) the choice of the plus sign is preferred when 
K > 0, and, conversely, for k < 0 it is necessary to choose the minus sign. 

Otherwise, in the most interesting region k 1 we receive a negative 
value for the 
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Thus 






+ 


Then 


Ur = 


2(l + /32) I + /32 

\/l + 






1 + /3' 


/S 


2 _ ^2 




(2.14.1.46) 


(2.14.1.3a) 


Obviously, that a small excess of the |k| over one means that Ur value will 
be small. 

|K| = l + e —> Ur = O{e){e>0) (2.14.1.9) 

Mgr Mgr 


A = 


We also have 


2KUr A 
-2 


= 0(e) > 0 


r = 1/91+0(e) 


(2.14.1.6a) 

(2.14.1.10) 


2.14.2 Sufficient conditions of stability for levitation 
of the Orbitron 


The expressions of A, B, C take the form 


4 _ ]\/r, 2 3 Mr„-Ij^u _ z JD _ 

B = —jiu'' Bz^r 
C = -nu^Bz,z 


‘^ 1 ^ 4 +Mrl (i'^^+^ 2)+M5(1-U) j 




Mg (e 


2 _ Mgr \ . 

4k2a1 ’ 


TTZr 


(2.14.2.1) 

Given (2.14.1.6a,2.14.1.9,2.14.1.10), we get sufficient conditions in the 
form 


'“ = Bi-4 = + B,,„) - 0(£) > 0; 

'>=l?;S = OWi 

<= = - 0 (f) > 0 ; 

ac — h^> 0; 

^ A = (jgiBz + cuTTo - /±a;^ > 


(2.14.2.2) 
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at that 0(e) > 0. 

From (2.14.2.2) and 0(e) > 0 it follows that geometric conditions are 
fulhlled 


^ 0 , 

{^Bz^r + Bz^rr) = 


-a 


-a 


-Bz^zz + f-Sz,r) > 0; 


(2.14.2.3) 


that ensure compliance with sufficient stability conditions (2.14.2.2) for suf- 
hciently small e > 0. 

Thus condition (2.14.2.3), together with dynamic condition 

ojttq S> —afiBz: + + Mgr (2.14.2.4) 


gives a full set of conditions of the system stability, i.e. Orbitron can levitates. 

Remark 7. Parameter to a certain extent is adjustable, as contains 
the contribution of the held B^ that does not affect on the equilibrium of 
the magnetic force and gravity force, and this parameter occures only in the 
(2.14.2.4). 
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